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Abstract — A pair of basis functions is presented for the 
surface integral, method of moment solution of scattering by 
arbitrarily-shaped, three-dimensional dielectric bodies. 
Equivalent surface currents are represented by orthogonal 
unit pulse vectors in conjunction with triangular patch 
modeling. The electric field integral equation is employed 
with closed geometries for dielectric bodies; the method may 
also be applied to conductors. Radar cross section results are 
shown for dielectric bodies having canonical spherical, 
cylindrical, and cubic shapes. Pulse basis function results are 
compared to results by other methods. 

Index Terms — electromagnetic scattering, basis 
functions, method of moments, dielectric, conducting, EFIE. 

I. Introduction 

This work extends the use of a pair of orthogonal basis 
functions introduced in [1] for the determination of 
surface currents on a three-dimensional, arbitrarily-shaped 
electromagnetic scatterer. In [1], a solution method was 
demonstrated for perfect electric conductors (PEC) using 
unit pulse basis vectors; in this work, the method is further 
developed to include dielectric bodies. 

Rao, Wilton, and Glisson (RWG) developed the RWG 
basis functions for method of moments (MoM) 
conducting body solutions in conjunction with triangular 
patch surface modeling [2]. Numerous other basis 
functions have been used for MoM solutions, including 
rooftop [3], wavelet [4], and Trintinalia-Ling [5] functions 
for specialized applications. However, for “all-purpose” 
modeling of arbitrarily-shaped, three-dimensional 
scatterers, the RWG basis functions remain in very wide 
usage. We propose that for lossy dielectric or dielectric 
(material) bodies, a pair of orthogonal basis functions 
should be more suitable in that the basis functions 
represent equivalent surface electric and magnetic currents 


also defined as orthogonal to each other in each triangular 
patch. Improved solution reliability and accuracy should 
result from orthogonal basis functions, which produce 
strongly diagonal, well conditioned impedance matrices. 
In addition, the pulse basis vector functions are 
particularly simple in concept and in definition, which 
makes them attractive to explore. Orthogonal pulse basis 
functions are easier to implement in code than, for 
example, it X RWG basis functions [6], ft defined as the 
unit outward normal vector. The paired pulse basis vector 
solution method may be extended to include collections of 
scatterers, composite dielectric/conducting bodies, or 
bodies composed of numerous regions of constant 
permittivity and permeability. 

II. Problem Formulation 
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Fig. 1. Arbitrarily-shaped, three-dimensional scatterer. 

Considering the dielectric scatterer to be illuminated by 
an external source, we label the outer and inner regions of 
the scatterer as 1 and 2, respectively, shown in Fig. 1. We 


apply the Equivalence Principle and enforce the continuity III. NUMERICAL SOLUTION PROCEDURE 

of the tangential E-field across the boundary surface S to 

write A. Basis Functions 
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just outside the surface. The superscript s and i refer to 
scattered and incident fields. J and M are the equivalent 
electric and magnetic surface currents residing on S. The 
scattered tangential electric field in (1) or (2) is given by 
the electric field integral equation (EFIE) as: 
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where A is the magnetic vector potential. 


A = -^~ f JGdS, (4) 
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0 is the electric scalar potential, 0 — V • A, (5) 
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F is the electric vector potential, 

F=—[ M-GdS, (6) 

An is 

e -3P R 

G is the Green’s function, G — , (7) 


Within each pair of adjacent triangular patches T„ + and 
T~, the basis vectors for J n and M„ are f, and g„, defined 
to be normal and parallel to the n ,h edge, respectively, as 
shown in Fig. 2. The basis vector functions have unit 
value and exist only in the shaded source regions, which 
need not be coplanar. These regions comprise two thirds 
of the area of the T„ + and T~ triangles. Plus and minus 
names have been assigned to designate the direction of 
electric current flow from plus to minus. 

n h edge 



the surface currents J„ and M„ will be calculated normal and 
parallel to the edge, respectively. 


At the n ,h edge. 


P is the phase constant 



( 8 ) 


J n= a n(fn +fn ) (9) 


R is the distance between source and observation points, s 
is electric permittivity, and ju is magnetic permeability. 

(1) and (2) are rewritten by expanding J and M in all 
vector and scalar potential terms as pulse basis functions 
and incorporating the regions’ respective values for s and 
/j . For a closed triangular mesh containing N edges, the 
expanded MoM matrix equation contains a 2 N x 2 N 
impedance matrix to solve for N values of J and N values 
of M. These surface currents may then be used to 
calculate the radar cross section (RCS) from any 
observation angle. 


M n=Yn(gn ) ( 10 ) 

where a n and y„ are the expansion coefficients to be 
determined. For the mn h impedance matrix term, vector 
potentials A mn and F m „ are evaluated by integrations (4) 
and (6) over the shaded current source regions. 

B. Testing Functions 

Testing vector functions at the m' h edge are defined as 
£,„ + and l m ~ , as shown in Fig. 3. £,„ + extends from the T m + 
centroid to the midpoint of the m"‘ edge; l m extends from 
the midpoint of the m‘ h edge to the T m ~ centroid. 
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Fig. 4. Observation points are the S m centroids for evaluation 
of A mn and F mn , the T m centroids for evaluation of @„ m . 


C. Testing Equations 

In each of regions 1 and 2 (shown in Fig. 1), N testing 
equations (11) are written in matrix form to solve (1) and 
(2) for J and M at the patch edges m = 1 to N. For region 
1, the right hand side of (11) has assigned nonzero values 
for the incident field; in region 2, the right hand side is 
zero because there is no incident field inside the scatterer. 
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[2] and [6] describe more details of triangular patch 
modeling with RWG basis functions; [1] further explains 
the implementation of unit pulse vector basis functions. 

D. Observation and Source Points for Potential 
Evaluation 

For the mn ,h impedance matrix term, vector potentials 
A mn and F mn are evaluated as observed from the centroids 
of the 8„f and S„f triangles. These triangles are defined 
by nodes at the T m and T m centroids and the endpoints of 
the m ,h edge, shown in Fig. 4. 

For the mn h impedance matrix term, scalar potential 
0 mn is evaluated as observed from the centroids of the T m + 
and T m triangles, shown in Fig. 4. The source region is 
the entirety of the T„ and T n ~ triangles shown in Fig. 2. 

E. Calculation of Potential Terms 

A is calculated by a straightforward application of (4) 
with G computed numerically. 

0 is calculated from (5), making use of the Divergence 
Theorem. We assume that, everywhere in T n , 0 n is the 
divergence of A through the three sides of the T n triangle. 


That divergence is calculated as 

V • A - (JjLj + J 2 L 2 + J 2 L 3 )lr GdS (12) 

where L ; , L 2 , and L } are the lengths of the sides of the T„ 
triangle and J h J 2 , and J 2 are the scalar current 
components normal to those sides. Application of the 
testing functions has eliminated the necessity to compute 

V 0 because the dot product of the testing vector with 

V 0 reduces to A Clover the testing path. 

V x — F is evaluated as follows: 
s 


Vxif = Vx — f M-GdS 
s 4/r J s 

= — f VGxMdS (13) 
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where S is the shaded area depicted in Fig. 2. This 
“deleted integral” formulation follows from evaluating the 
Cauchy principal value for the singular self terms when, 
as in our case of 8 m centroid evaluation, the observation 
point is not on a patch edge [7]. 

F. Incident E-field 

The incident E-field is evaluated at the S„, + and 
8 m ~ centroids, which are located two thirds of the way 
along the t n f testing path and one third of the way along 
the testing path. 
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IV. Numerical Results 

Examples are shown of scattering calculations for the 
canonical shapes of a sphere, a finite cylinder, and a cube. 
These geometries are represented in Fig. 5. In each case, 
an incident plane wave is assumed to be x-polarized, 
traveling in the +z direction. The pulse basis vector 
results are compared to results of calculations by other 
methods. 





Fig. 5. RCS has been calculated for a dielectric sphere, a 
dielectric cylinder, and a dielectric cube separately. 

A. Sphere 

A dielectric sphere having diameter = 0.2 X and £r = 4 
is meshed with 500 triangles (750 edges) as shown in Fig. 
5. The bistatic RCS is shown versus 6 for observation 
angles ($=0 (Fig. 6) and ($=90 (Fig. 7). These results 
are compared with the analytic Mie series results [8], 

B. Cylinder 

A dielectric cylinder having diameter = 0.2 X, height = 
0.2 X, and s R = 3 is meshed with 320 triangles (480 edges) 
as shown in Fig. 5. The bistatic RCS is shown versus 6 
in Figs. 8 and 9. These results are compared with body of 
revolution (BOR) MoM results. 

C. Cube 

A dielectric cube having side length = 0.2 X and e R = 4 
is meshed with 480 triangles (720 edges) as shown in 
Fig. 5. The bistatic RCS is shown versus 0 in Figs. 10 
and 11. These results are compared with surface integral 
equation (SIE) and volume integral equation (VIE) MoM 
results reported by Sarkar et al. [9], Sarkar’s SIE 
formulation used 192 triangles (288 edges), while the VIE 
formulation used 512 cubes. 


D. Discussion 

For the examples of dielectric sphere, cylinder, and 
cube, the paired pulse basis expansion functions produce 
RCS values very comparable to those values determined 
by other accepted methods. As the modeling mesh is 
made finer to increase spatial resolution, pulse results 
converge to the expected values. Although real values of 
s have been used in the examples, complex values may be 
used to model a lossy dielectric or conductive material. 
The paired pulse basis vector solution method may also be 
extended to include collections of scatterers, composite 
dielectric/conducting bodies, or bodies composed of 
numerous regions of constant permittivity and 
permeability. 

V. Conclusion 

A pair of orthogonal unit pulse basis vector functions 
has been applied to the solution of MoM scattering 
problems for three-dimensional dielectric objects in the 
shapes of a sphere, a finite cylinder, and a cube. Good 
agreement was obtained between these results and the 
results from other analytic and numerical methods. These 
results indicate that the paired pulse basis vector method 
can be used for arbitrarily-shaped bodies, given a 
triangular mesh resolution appropriate for the electrical 
size of the body. The orthogonal spatial relationship 
between the electric and magnetic current basis functions 
should provide stable solutions for problems involving 
single, multiple, or hybrid dielectric/conducting bodies. 
Although other methods may give more accurate results 
for specific geometries, such as the Mie series method for 
a sphere or the BOR method for a cylinder, the paired 
pulse basis function method in combination with 
triangular surface patches will give good answers for 
many different geometries with fewer computations than 
the volume integration method. 
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Fig. 6. Bistatic RCS at (j>= 0° for a dielectric sphere, S R = 4. 
Pulse method (500 triangles) vs. Mie series method. 



Fig. 8. Bistatic RCS at (j>= 0° for a dielectric cylinder, Sr = 3. 
Pulse method (320 triangles) vs. BOR method. 



Fig. 10. Bistatic RCS at </> = 0° for a dielectric cube, S R = 4. 
Pulse method (480 triangles) vs. Sarkar’s SIE method (192 
triangles) and Sarkar’s VIE method (512 cubes) [8], 



co 

o 


0 30 60 90 120 150 180 

Theta [degrees] 

Fig. 7. Bistatic RCS at (j>= 90° for a dielectric sphere, £r = 4. 
Pulse method (500 triangles) vs. Mie series method. 



Fig. 9. Bistatic RCS at <j>= 90° for a dielectric cylinder, s R = 
3. Pulse method (320 triangles) vs. BOR method. 
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Fig. 11. Bistatic RCS at </)= 90° for a dielectric cube, e R = 4. 
Pulse method (480 triangles) vs. Sarkar’s SIE method (192 
triangles) and Sarkar’s VIE method (512 cubes) [8], 
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